Critical spectrum of fluctuations for deconfinement at proto-neutron star cores 
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We study the deconfinement of hadronic matter into quark matter in a protoneutron star focus- 
ing on the effects of the finite size on the formation of just-deconfined color superconducting quark 
droplets embedded in the hadronic environment. The hadronic phase is modeled by the non-linear 
Walecka model at finite temperature including the baryon octet and neutrino trapping. For quark 
matter we use an SU (3) / Nambu- Jona-Lasinio model including color superconductivity. The finite 
size effects on the just deconfined droplets are considered in the frame of the multiple reflection 
expansion. In addition, we consider that just deconfined quark matter is transitorily out of equi- 
librium respect to weak interaction, and we impose color neutrality and flavor conservation during 
the transition. We calculate self-consistently the surface tension and curvature energy density of 
the quark hadron inter-phase and find that it is larger than the values typically assumed in the 
literature. The transition density is calculated for drops of different sizes, and at different temper- 
atures and neutrino trapping conditions. Then, we show that energy-density fluctuations are much 
more relevant for deconfinement than temperature and neutrino density fluctuations. We calculate 
the critical size spectrum of energy-density fluctuations that allows deconfinement as well as the 
nucleation rate of each critical bubble. We find that drops with any radii smaller than 800 fm can 
be formed at a huge rate when matter achieves the bulk transition limit of 5—6 times the nuclear 
saturation density. 

PACS numbers: 12.39.Fe, 25.75.Nq, 26.60.Kp 



I. INTRODUCTION 

The cores of massive neutron stars have densities that 
may favor the nucleation of small droplets of quark mat- 
ter, that under appropriate conditions may grow con- 
verting a large part of the star into quark matter (for 
recent work see [IHlj and references therein). Since the 
the equation of state above the nuclear saturation density 
is still uncertain the problem is usually analyzed within 
a two-phase description in which a hadronic model valid 
around the nuclear saturation density po is extrapolated 
to larger densities and a quark model that is expected 
to be valid only for asymptotically large densities is ex- 
trapolated downwards. Within this kind of analysis some 
work has been performed recently in order to determine 
the effect of different hadronic and quark equations of 
state, as well as the effect of temperature, neutrino trap- 
ping and color superconductivity 

An important characteristic of the deconfinement tran- 
sition in neutron stars, is that quark and lepton flavors 
must be conserved during the process leaving just decon- 
fined quark matter that is transitorily out of equilibrium 
with respect to weak interactions (see [l[ and references 
therein). When color superconductivity is included in 
the analysis together with flavor conservation, it is found 
that the most likely configuration of the just deconfined 
phase is two-flavor color superconductor (2SC) provided 
the pairing gap is large enough Q. 

Since our previous studies were made in bulk, we shall 
now include finite size effects in the description of the just 



deconfined drops. To this end we employ the Nambu- 
Jona-Lasinio (NJL) model for quark matter and include 
finite size effects within the multiple reflection expan- 
sion (MRE) framework (sec [6 8] and references therein). 
Through this analysis we can determine the density of 
hadronic matter at which deconfinement is possible for 
different radii of the just formed quark drops using typ- 
ical conditions expected in the interior of protoneutron 
stars (PNSs), i.e. temperatures in the range of 0—60 MeV 
and chemical potentials of the trapped neutrino gas up 
to 200 MeV. Notice that to form a finite size drop it is 
needed some over-density with respect to the bulk tran- 
sition density in order to compensate the surface and 
curvature energy cost. Since this energy cost depends 
on the drop radius R, so does the necessary over-density 
necessary to nucleate it. Thus, we can derive a critical 
fluctuation spectrum (5p/p versus R) delimiting which 
over-densities can deconfine and grow unlimitedly and 
which ones will shrink back to hadronic matter. 

The paper is organized as follows. In Section II we 
briefly describe main features of the hadronic model. In 
Sect. Ill we present the quark model we use taking into 
account finite size effects and derive self-consistently the 
surface tension and the curvature energy from the MRE 
thermodynamic potential. In Sect IV we the find the den- 
sity at which color superconducting droplets of different 
radii can be nucleated in protoneutron star matter. In 
Section V we find the critical spectrum of density fluc- 
tuations as well as the nucleation rate of the critical size 
droplets. In Sect. V we present our conclusions. 
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II. THE HADRONIC PHASE 

The hadronic phase is modeled by the non-linear 
Walecka model (NLWM) [iO] including the whole 
baryon octet. The Lagrangian of the model is given by 



C = C, 



C 



M 



(1) 



where the indices B, M and L refer to baryons, mesons 
and lcptons respectively. For the baryons we have 



B 



0i 



7 M (idn - 9u>b u p - g pB f ■ p^) 



-(m B - g a B a) 



ipB, (2) 



with B extending over the nucleons N = n, p and the 
following hypcrons H = A, £ + , £°, E~, S~, and S°. 
The contribution of the mesons a, ui and p is given by 

1 be 



-i w MV + i to 2 , ujf, uj^ 



(3) 



where the coupling constants are 

9crB = XaB 9cr, 9luB = %uB 9u , 9pB = X pB 9p (4) 

with x g b, %uB and x p b equal to 1 for the nucleons and 
0.9 for hyperons. Electrons and neutrinos are included 
as a free Fermi gas in chemical equilibrium with all other 
particles. 

There are five constants in the model that are deter- 
mined by the properties of nuclear matter, three that de- 
termine the nuclcon couplings to the scalar, vector and 
vector-isovector mesons g (7 /'m (7 , g^/m^, g p /m p , and two 
that determine the scalar self interactions b and c. It 
is assumed that all hyperons in the octet have the same 
coupling than the A. These couplings are expressed as a 
ratio to the nucleon couplings mentioned above, that we 
thus simply denote x a , x u and x p . 

In the present work we use the following values for 
the constants: (g a /m a ) 2 = 11.79 fm 2 , (</ w /to w ) 2 = 



0.002947 and 



7.149 fm 2 , (g p /m p ) 2 = 4.411 fm 2 , b 
c = -0.001070. This model has been labelled as GM4 
in previous work [l[ . The equation of state is rather stiff 
and gives a maximum mass of 2 M Q , that seems to be 
more adequate in light of the recently determined mass 
of the pulsar PSR J1614-2230 with M = 1.97 ± 0.04M o 
[r|. This is the largest mass reported ever for a pulsar 
with a high precision. For more details on the explicit 
form of the equation of state the reader is referred to 
Ref. Q. 



the inclusion quark-quark interactions, which are the re- 
sponsible for color superconductivity. 

The corresponding Lagrangian is given by 



C = ijj (ift — m) ip 

8 



a=0 



(5) 



+2H ^2 [{i> ilhTA^A' ipe) (i>c ilbT A \ A , if)] 

A,A'=2,5,7 

where m = diag(m u , rrid, m s ) is the current mass ma- 
trix in flavor space. In the present paper we work in the 
isospin symmetric limit m u = ma = to. Moreover, n 
and Ai with i = 1, .., 8 are the Gell-Mann matrices corre- 
sponding to the flavor and color groups respectively, and 
To = \/2/3 If. Finally, the charge conjugate spinors are 
defined as follows: ipe = C 4> T and ipe = "0 T C: where 
%p = tp'j is the Dirac conjugate spinor and C = i"f 2 ~f°. 

In order to determine the relevant thermodynamic 
quantities we have to obtain the grand canonical thermo- 
dynamic potential at finite temperature T and chemical 
potentials pf c . Here, / = (u,d,s) and c = (r, g, b) de- 
notes flavor and color indices respectively. In the follow- 
ing we present the thermodynamic potential for the bulk 
system and then, in Sect. Ill B, we discuss the effects 
of finite size in the effective potential, for the spherical 
droplets. 



A. Quark matter in bulk 

Starting from Eq. ([5|), we perform the usual bosoniza- 
tion of the theory. We introduce the scalar and pseu- 
doscalar meson fields a a and ir a respectively, together 
with the bosonic diquark field A a- In what follows we 
will work within the mean field approximation (MFA), 
in which these bosonic fields arc expanded around their 
vacuum expectation values and the corresponding fluc- 
tuations neglected. For the meson fields this implies 
v = o a T a = diag(cr„,cr d ,er s ) and 7r a = 0. Concern- 
ing the diquark mean field, we will assume that in the 
density region of interest only the 2SC phase might be 
relevant. Thus, we adopt the ansatz Ag = A7 = 0, 
A2 = A. Integrating out the quark fields and working 
in the framework of the Matsubara and Nambu-Gorkov 
formalism wc obtain the following MFA quark thermo- 
dynamic potential fl^ IFA (T, p,f c , a u , ad, cr s , |A|) per unit 
volume (further calculation details can be found in Refs. 



III. THE QUARK MATTER PHASE 

The just deconfined quark matter phase is imple- 
mented by using the SU(3)f NJL effective model with 



r\MFA 

3 

V 



k 2 dk x - 

i=l 



u{xi,yi) 



AG 



(a 2 + a 2 + a 2 ) + 



2H 



(6) 
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where A is the cut-off of the model and u>(x, y) is defined 
by 



u{x,y) 



4-Tln[l 
-TIn[l -+ 



-(?-v)/T-\ 



-(x 



(7) 



with 



2=1,2 



2=6,7 



3=8,9 

Vi 

Vi 

2/6,7 



3=3,4,5 




(8) 



Here, E = \Jk 2 + M 2 and E s = ^Jk 2 + M|, where Mj = 
rrif + <jf. Note that in the isospin limit we are working 
o'u = 0"d = a and, thus, M M = = M. 

The total thermodynamic potential of the quark mat- 
ter phase (super-index Q) is obtained by adding to {l^ IFA 
the contribution of the leptons and a vacuum constant. 
Namely, 



= n A J FA + ^ 



(9) 



where f2 e and are the thermodynamic potentials of 
the electrons and neutrinos respectively. For them we 
use the expression corresponding to a free gas of ultra- 
relativistic fermions 



V 



-Pi 



P± 
24tt 2 



tfT 2 



7tt 2 T 4 



12 360 

" (10) 

where P stands for the pressure, I = e,v ei and the de- 
generacy factor is j e = 2 for electrons and 7^ = 1 for 
neutrinos. Notice that in Eq. © we have subtracted 
the constant f2 vac = — P mc V in order to have a vanishing 
pressure at vanishing temperature and chemical poten- 
tials. More details are shown below in Sec. IIII CI 



B. Finite size effects: inclusion of MRE 



where 



and 



fs — 

rrif 



1/2 k 

1 arctan 

8n \ 7r rrif 



— (12) 



fa — = 
rrif 



12tt 2 



3k f n k 

1 arctan 

2m f \ 2 rrif 



(13) 



are the surface and curvature contributions to the new 
density of states respectively. For fc we employ the Mad- 
sen ansatz Q because its functional form for any finite 
quark mass has not yet been derived in the MRE frame. 

As shown in [H , the density of states of MRE for mas- 
sive quarks is reduced compared with the bulk one, and 
for a range of small momentum becomes negative. The 
way of excluding this non physical negative density of 
states is to introduce an infrared cutoff in momentum 
space (see @ for details). Thus, in our quark model, 
which includes the MRE formalism for finite spherical 
droplets of color superconducting quark matter, we have 
to perform the following replacement 



k 2 dk 
2ir 2 



k 2 dk 



■PMRE- 



(14) 



In order to obtain the value of K-ir, we have to solve 
the equation pMRE = with respect to the momentum 
k and take the larger root as the IR cut-off. Note that 
Pmre depends on the quark mass and the radius of the 
droplets and consequently the Am has the same depen- 
dence (more details are given below in Sec. IIII Cp . 
Therefore, the full thermodynamic potential reads: 



AIRE 

V 



/ k 2 dk A , 

i 1 tjl i J2 i JZ\ | l^| 2 

+^K + a d + a s ) + — 



(15) 



Multiplying on both sides of the last equation by the 
volume of the quark matter drop and rearranging terms 
we arrive to the following form for Q MRE 



Sl Q MRE = -P Q V + aS + 1 C, 
where the pressure P® is given by 



(16) 



In the present work we consider the formation of fi- 
nite size droplets of quark matter. The effect of finite 
size is included in the thermodynamic potential adopting 
the formalism of multiple reflection expansion (MRE; see 
Refs. [||-[8| and references therein). 

In the MRE framework, the modified density of states 
of a finite spherical droplet is given by 



PMRE(k,m f ,R) = 1 



— fs (— 

kR \1Tlf 



127T 2 

(fci?) : 



-fc 



k 

rrif 
(11 



pQ — _ uli MRE 



8V 



the surface tension is 



on 



MRE 



dS 



^{x uVl ) 

T,n,S,C J ^iR z " i-i 

1/2, J.J^ l A | 2 



T,fi,V,C 



+Pe + P„ e ~ Pvac 



/•A 9 
2 / kdkf s y^ y uj(xi, 
J^ir i= i 



(17) 



Vi), (18) 
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m u , d [MeV] 


m s [MeV] A [MeV] 


GA 2 


H/G 


Set 1 


5.5 


112.0 602.3 


4.638 


3/4 


Set 2 


5.5 


110.05 631.4 


4.370 


3/4 



TABLE I: The two sets of NJL parameters, 
and the curvature energy density is 

rA 9 



1 



AIRE 



dC 



T,n,V,S 



dkfc^2^(xi,yi). (19) 

»=i 



Here we are considering a spherical drop, i.e. the area is 
S = AttR 2 and the curvature is C — 8irR. 

From the grand thermodynamic potential &mre we 
can readily obtain the number density of quarks of 
each flavor and color n/ c = —V~ 1 dQ® IRE /d[if c , the 
number density of electrons n e = —V^ 1 dfl < ^ IRE /dfi e , 
and the number density of electron neutrinos n Ve = 
— V~ 1 dfl C j' [RE /dfiv a . The corresponding number densi- 
ties of each flavor, n/, and of each color, n c , in the 
quark phase are given by nf =J2c n fc an< ^ n c = J2f n fc 
respectively. The baryon number density reads tib = 
lE/c n /c = ( n « n d + n s )/3. Finally, the Gibbs free 
energy per baryon is 



^~ [ X! ^f c n f c + t le Ue + Vv e n Ve | . (20) 



C. Parametrizations 

For the NJL model, the values of the quark masses and 
the coupling constant G can be obtained from the meson 
properties in the vacuum. Here we use a set of parameters 
taken from [l8|, but without the 't Hooft flavor mixing 
interaction. The procedure, obtained from [lj| , is to keep 
A and m fixed and then tune the remaining parameters G 
and m s in order to reproduce M — 367.6 MeV and M s = 
549.5 MeV at zero temperature and density. An estimate 
of H/G can be obtained from Ficrz transformation of the 
one-gluon-exchange interactions in which case one gets 
H/G = 0.75. The resulting parameter sets are given in 
Table fl] For this set of parameters we get f2 vac /F = 
-P vac = -4301 MeV/fm 3 and -5099 MeV/fm 3 (for set 
1 and 2 respectively). 

As we previously mentioned, the value of Aju is the 
largest root when solving pmre = with respect to k, 
depending on mt and R. For a given mass, we find nu- 
merically that Air(R) can be fitted as follows 



Air = aR b 



(21) 



with R in fm and Ajr in MeV. The coefficients a and 
b are found to be 135.45 and —0.85 for m u = 5.5 MeV; 
227.72 and -0.86 for m s = 110.5 MeV; 228.60 and -0.87 
for m s = 112 MeV. 



D. Flavor conservation, color neutrality and other 
conditions 

In order to derive a quark matter EOS from the above 
formulae it is necessary to impose a suitable number of 
conditions on the variables {/U/ c }, /i e , Hv c , c, a s and A. 
Three of these conditions are consequences from the fact 
that the thermodynamically consistent solutions corre- 
spond to the stationary points of ^i^f RE with respect to 
a, (7 S , and A. Thus, we have 



on 



AIRE 



= 0, 



AIRE 



da. 



0, 



u * ''AIRE 

d\A\ 



0. (22) 



To obtain the remaining conditions one must specify the 
physical situation in which one is interested in. As in 
previous works P, H, Hlj , we are dealing here with 
just deconfined quark matter that is temporarily out of 
chemical equilibrium under weak interactions. The ap- 
propriate condition in this case is flavor conservation be- 
tween hadronic and deconfined quark matter. This can 
be written as 



\rR \rQ 



f = u, d, s, e, v e 



(23) 



being Yf = nf/n 1 ^ and Y i 



Q 



Q /„Q 



i f /n 



the abundances 



of each particle in the hadron and quark phase respec- 
tively. It means that the just deconfined quark phase 
must have the same "flavor" composition than the /3- 
stable hadronic phase from which it has been originated. 
Notice that, since the hadronic phase is assumed to be 
electrically neutral, flavor conservation ensures automat- 
ically the charge neutrality of the just deconfined quark 
phase. The conditions given in Eq. (|23|) can be combined 
to obtain 



rid 

lis 

3n» 



V n u , 

2n u — nd — n s , 



(24) 
(25) 
(26) 
(27) 



where n.; is the particle number density of the i-species 
in the quark phase. The quantities £ = Yf /Yf , r\ = 
Yf /Yf and k = Yf /Yf are functions of the pressure 
and temperature, and they characterize the composition 
of the hadronic phase. These expressions are valid for 
any hadronic EOS. For hadronic matter containing n, p, 
A. E+, S°, E- 3- and 5°, we have 



n p + 2n n + n A + n s o + 2n s - + n s - 
2n p + n n + riA + 2n E + + n s o + n s o ' 
n\ + n s + + n E c + n s - + 2n s o + 2n s - 
2n p + n n + Ua + 2n s + + n s o + n s o 

< 

2n p + n n + n\ + 2n s + + n s o + nu° 



(28) 
(29) 



(30) 



Additionally, the deconfined phase must be locally col- 
orless; thus it must be composed by an equal number of 
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red, green and blue quarks 

TL r = 11 g = ni,. 



(31) 



Also, ur, ug, dr, and dg pairing will happen provided 
that |A| is nonzero, leading to 



nd a , 



II 



ndr- 



(32) 



In order to have all Fermi levels at the same value, we 
consider Q 



n 



n s b 



I I .s 



(33) 



These two equations, together with Eqs. ([31]) and (|32|) 
imply that n ur = n ug = n dr = n dg and n sr = n sg = n sb 

a. 

Finally, including the conditions Eqs. (|22p we have 13 
equations involving the 14 unknowns (c, a Sl |A|, /i e , p, Ve 
and {/i/ c })- For given value of one of the chemical poten- 
tials (e.g. /J, ur ), the set of equations can be solved once 
the values of the parameters £, ry, k, the temperature T 
and the radius of the drop R are given. Instead of p, ur , we 
can provide a value of the Gibbs free energy per baryon 
SVark and solve simultaneously Eqs. ([24"f - (f33"|) together 
with Eqs. (|22p in order to obtain cr, <r si |A|, p e , fi Ve and 

{M/c}- 



IV. DECONFINEMENT OF COLOR 
SUPERCONDUCTING DROPLETS IN 
PROTO-NEUTRON STAR MATTER 

The total thermodynamic potential of a quark matter 
drop immersed in an homogeneous hadronic environment 
is Q = fl H + n® IRE , where the indexes H and Q refer 
to the hadronic and the quark phase respectively. For 
the hadronic phase we have n H = —P H V H and for the 
quark phase we have n^ IRE = —P^V® + aS + 7C. 

The condition of mechanical equilibrium is given by 



an 

dR 



dn 



T,n,V 



OR 



dn 



MRE 



T,n,V 



dR 



= 0. 



T,n,V 



Since V = V® + V H is constant, we have: 



n ndVQ dS dC nH dV Q 

P w h a h y h P 

dR dR ' dR dR 



0. 



(34) 



(35) 



Thus, for a spherical droplet the condition for mechan- 
ical equilibrium reads (c.f. Q): 



pQ _ 2a _ 27 _ H 
R R 2 



0. 



(36) 



Notice that in the bulk limit R — > 00 we find the standard 
condition P H = P Q . 

Additionally, we assume thermal and chemical equilib- 
rium, i.e. the Gibbs free energy per baryon are the same 



for both hadronic matter and quark matter at a given 
common temperature. Thus, we have 



fJ 



T 



H 



(37) 



If we fix the radius R of the deconfined drop for a given 
temperature T H and neutrino chemical potential of the 
trapped neutrinos in the hadronic phase fj,^ , there is an 
unique hadronic pressure P H at which the equilibrium 
conditions are fulfilled. Notice that, differently than in 
the bulk case studied in pj, P H and P® arc not equal 
because of the appearance of a surface and a curvature 
term in the condition for mechanical equilibrium. Also, 
the mass-energy density p H and at the equilibrium 
point are different in general. Similarly, while the abun- 
dance Y Ve of neutrinos is the same in both the hadronic 
and just deconfined quark phases, the chemical potentials 
/xj and fj,^ are different. 

In the present work we want to describe thermody- 
namic conditions analogous to those encountered in pro- 
toncutron stars. From numerical simulations of the first 
tens of seconds of evolution [23j], we know that the pro- 
toncutron star cools from T ~ 40 MeV to temperatures 
below 2-4 MeV in about one minute. In that period, 
the chemical potential ii„ of the trapped neutrinos de- 
creases from ~ 200 MeV to essentially zero. Then, tak- 
ing into account typical protoneutron star conditions, we 
have solved Eqs. (|3"6"|) and (|3"T|) together with the condi- 
tions showed in Sections II and III, for temperatures in 
the range 0-60 MeV, (1% fii the range - 200 MeV and 
for different radii of the color superconducting droplets. 
The results are displayed in Figs. 1 and 2 for the param- 
eterizations of the equations of state given in previous 
sections. 

In Fig. 1 we show the mass-energy density of hadronic 
matter (in units of the nuclear saturation density po = 
2.7 x 10 g cm 3 ) above which it is energetically favorable 
the conversion of a portion of hadronic matter into an 
identical deconfined drop (i.e. with the same T, P, g, 
Yf and radius R) . For the calculations we considered the 
following possibilities: set 1 of the NJL model in the left 
panels and set 2 in the right panels, temperatures in the 
range 0-60 MeV, and neutrino chemical potential p^f = 
0, 100, 200 MeV. Within each plot we consider different 
values of R; from right to left the curves correspond to 
R[fm] = 2, 5, 10, 20, 50, 500, 00. Notice that, in the six 
panels of Fig. 1, the curves for the bulk case (R = 00) 
are almost coincident with the curves for R = 500 fm. 
The main observed effect in Fig. 1 is that the transition 
density increases considerably for small radii. This is 
a consequence of the significant increase of the surface 
tension a and the curvature energy 7 for small R (sec 
Fig. 2). 

In Fig. 2 we display the surface tension a and the cur- 
vature energy 7 as a function of temperature for different 
radii of the droplet and with the set of {/i/ c , <r, (Tf, A} 
that arises from the equilibrium condition discussed be- 
fore. Notice that a and 7 do not depend explicitly on the 
radius R of the drop (see Eqs. (jT5)) and (fH)]0 . but for a 
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NJL setl NJL set2 




FIG. 1: For different temperatures we show the mass-energy density p of hadronic matter above which it is energetically 
favorable the conversion of the hadronic phase into deconfined drops of radius R. po — 2.7 x 10 14 g cm 3 is the nuclear saturation 
density. For the hadronic phase we employ the model introduced in Sec. II. For the quark phase we use set 1 of the NJL 
model presented in Sec III in the left panels, and set 2 in the right panels. In each figure the neutrino chemical potential in the 
hadronic phase is fixed to the following values: = 0, 100, 200 MeV. Within each plot, each curve correspond to a different 
value of R. From right to left we have _R[fm] — 2, 5, 10, 20, 50, 500, oo. Notice that in all figures the curves for the bulk case 
(R — oo) are almost coincident with the curves for R = 500 fm. A possible interpretation of these figures is best illustrated 
with an example: let us consider a hadronic system with a density p H ~ 6 x po at a temperature T = 10 MeV and with 
trapped neutrinos having a chemical potential p Ve = 200 MeV. This state is represented by a circle in the two upper panels 
of this figure. Looking at the position of the circle with respect to the curves corresponding to different radii, we see that it 
is energetically favorable to convert hadronic matter into quark drops provided they have a radii R > 50 fm (if quark matter 
is described by the NJL set 1 EoS) and R > 10 fm (using NJL set 2 EoS). That is, radii whose curves lie to the left of the 
circle are energetically favored while those on the right are forbidden because of the energy cost of forming a surface. A more 
realistic interpretation of the figure is given in Sec. V. 



given R, T and p^ there is a unique set of {p,f c , cr, a/, A} 
satisfying the equilibrium conditions, i.e. a and 7 depend 
on R through {pf c , cr, er/, A}. From Fig. 2 we see that 
while both a and 7 arc roughly independent of the tem- 
perature, there is a larger dependence with the radius 
R leading to the behavior of the transition density ex- 
plained in the caption of Fig. 1. Notice also that the 
here-found values of a are larger that those found within 
the MIT Bag model which are typically ~ 30-60 MeV 
fm -2 for drops larger than few fm [24| . 



V. FLUCTUATIONS AND DECONFINEMENT 

According to the theory of homogeneous nucleation, 
the free energy involved in the formation of a spherical 
quark bubble of radius R is given by (25j 

AO = R 3 AP + AiraR 2 + 8irjR, (38) 

3 

where AP = P Q - P H is the pressure difference between 
internal and external parts of the bubble. For given AP, 
a and 7, the extremal points (maximum or minimum) of 
AO are obtained from dAil/dR = 0, which leads to Eq. 
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FIG. 2: The surface tension a and the curvature energy 7 as a function of temperature for different radii of the droplet and 
with the set of {/x/ c , cr, 07, A} that arises from the equilibrium conditions leading to Fig. 1. Different curves correspond to 
different radii R; in the case of a we have i?[fm] = 2, 5, 10, 20, 50, 500, 00 from bottom to top while in the case of 7 the values are 
the same but from top to bottom. Both a and 7 are roughly independent of the temperature but there is a larger dependence 
with R. 



(|36| . Thus, the critical radii are given by: 

R± = £p{ 1± ^ T + l ) . ( 39 ) 

with b = 2-fAP/a 2 . 

For b < — 1 both solutions are complex and Ail is a 
monotonically decreasing function of R. This means that 
any small fluctuation of one phase into the other will 
gain energy by expanding and a rapid phase transition 
is likely to occur. For b > — 1, AO has a local minimum 
at i?_ and a local maximum at R + . For Af2(i?+) < 
we have again that any small fluctuation is energetically 
favored. For AQ(i£+) > bubbles with radii larger than 
i?+ gain energy by growing unlimitcdly, while those be- 
low the critical size gain energy by shrinking to zero (if 
R- < 0) or to i?_ (if i?_ > ). In this case, the stan- 
dard assumption in the theory of bubble nucleation in 
first order phase transitions is that bubbles form with a 
critical radius R + [2Fjj . 

The approach wc adopted in the previous section is 
closely related to what we explained in the above para- 
graph. Instead of finding the critical radius for arbitrary 



values of AP, a and 7, we fixed R and found the corre- 
sponding AP, a and 7 that satisfy the conditions pre- 
sented in Sees. II, III and IV. Since Eq. (|36|) is satisfied 
by construction, the radius R is precisely the critical ra- 
dius R+ introduced in Eq. (|39[) . given that we choose the 
solution that verifies d 2 A£l/dR 2 < 0. 

In the light of the previous discussion we may give an- 
other interpretation to the results presented in Fig. 1. 
Let us consider, for simplicity, a uniform hadronic sys- 
tem with trapped neutrinos having a chemical potential 
fj,^ and a constant mass-energy density infinitesimally 
to the right of the curve with R = 00 in the T — p H /p 
plane corresponding to the same /i^. Since this density 
corresponds to the bulk transition we shall refer to it as 
Pbulk C^i l 1 ^ ) ■ F° r sucn a density, it is favorable to con- 
vert the system into quark matter in bulk, but this is 
not possible in practice due to the surface and curvature 
energy cost. 

However, fluctuations in the independent thermody- 
namic variables {T,p H ,^} of the hadronic fluid may 
drive some portion of it to a state described by {T + 
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FIG. 3: Critical spectrum for fluctuations that allow deconfinement of hadronic matter in a protoneutron star. As in the 
previous figures, the results were calculated using GM4 + NJL setl on the left panel and GM4 + NJL set 2 on the right panel. 
The spectrum does not depend significantly on the temperature and the chemical potential of trapped neutrinos. Fluctuations 
in hadronic matter having a given 8p H / p H are able to grow if they have a size R larger than the here-shown critical one. 



ST, p H + 5p H ,pH + SpHj. Nevertheless, notice that the 
curves of Fig. 1 are quite vertical, i.e. the transition is 
not very sensitive to changes in T. The same holds for 
variations in pjj . Thus, we shall consider only energy- 
density fluctuations with radius R fi that drive some part 
of the hadron fluid to a density p^ = p^ aik + Sp H to the 
right of a curve with a given R. As explained in the 
caption of Fig. 1, quark drops with radii larger than R 
are energetically favored. Thus, if the fluctuation has 
a size Rfi larger than R it will be energetically favor- 
able for it to grow indefinitely. In order to quantify this, 
we calculate the difference 5p H between p^ uik and the 
hadronic density of the point that allows nucleation with 
radius R or larger. In such a way we can construct a 
critical spectrum Sp H /p H as a function of R for differ- 
ent values of T and p^ as seen in Fig. 3. Fluctuations 
of a given over-density Sp H / p H must have a size larger 
than the critical value given in Fig. 3 in order to grow. 
Equivalently, fluctuations of a given size must have an 
over-density 5p H / p H larger than the critical one for that 
size. 

We can calculate the formation rate of critical bubbles 



through 

r « T 4 cxp{- Sn c /T). (40) 

where in our case <5f2 c is the work required to form a quark 
bubble with the critical radius from hadronic matter at 
the bulk transition point 

A'jr 

6Q C = -—R 3 (PQ - P b H ulk ) + 4naR 2 + 8^ 7 i?, (41) 

Instead of T 4 , different prefactors are used for T in other 
works (see e.g. [![). However, this fact does not affect 
significantly the results because T is largely dominated by 
the exponent in Eq. (|40|) ; i.e. we always have log 10 T w 
log 10 (prefactor) — Sfl c /[T ln(10)] with the second term 
much larger than the first. 

The results are given in Fig. 4 and show that critical 
bubbles with R > 800 fm are strongly disfavored while 
those with R < 800 fm have a huge rate. In practical 
situations, i.e. at neutron star cores, this means that 
if fluctuations lead hadronic matter to the bulk transi- 
tion point, quark drops with R < 800 fm will nucleate 
instantaneously. 
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FIG. 4: Nucleation rate for bubbles of the critical size. 



VI. SUMMARY AND CONCLUSIONS 

In the present paper we have studied the deconfinc- 
ment of quark matter in protoneutron stars employing a 
two-phase description of the first order phase transition. 
For the hadronic phase we used the non-linear Walecka 
model (Sect. II) which includes the whole baryon octet, 
electrons and trapped electron neutrinos in equilibrium 
under weak interactions. For the just deconfined quark 
matter we used a SU (3) f NJL model including color su- 
perconducting quark-quark interactions. In this phase, 
finite size effects are included via the multiple reflection 
expansion framework (Sect. III). Wc considered that 
color superconducting quark droplets are formed in me- 
chanical and thermal equilibrium with the hadronic en- 
vironment when the Gibbs free energy of both phases are 
equal. Also, since deconfmement is a strong interaction 
process, we consider that the just formed quark phase has 
the same flavor composition than the hadronic /3-stable 
phase, and consequently it is out of chemical equilibrium 
under weak interactions. The further /3-cquilibration of 
such quark phase is not addressed within the present pa- 
per. 

Through this model we determined the density of 
hadronic matter at which deconfincmcnt is possible for 
different radii of the just formed quark drops, as well as 



for different values of the temperature and the chemical 
potential of the trapped neutrinos (see Fig. 1). We have 
used typical conditions expected in the interior of pro- 
toneutron stars, i.e. temperatures in the range of — 60 
MeV and chemical potentials of the trapped neutrino gas 
up to 200 MeV. As expected, we find that the transi- 
tion density increases significantly when the radius of the 
droplet decreases. This is a consequence of the increase 
of the surface tension a and the curvature energy 7 for 
small radii. Additionally, we found that a and 7 at the 
deconfmement point are almost independent of the tem- 
perature and the neutrino chemical potential (see Fig. 
2). 

When the bulk transition density is reached in the core 
of a neutron star, it is energetically favored to convert 
a macroscopically large portion of hadronic matter into 
quark matter. But in practice, it is needed some over- 
density with respect to the bulk transition density in or- 
der to compensate the surface and curvature energy cost 
of a finite drop. Since this energy cost depends on the 
drop radius, so does the necessary over-density and over- 
pressure necessary to nucleate it. Thus, we can derive a 
critical fluctuation spectrum Sp H / p H versus R delimit- 
ing which fluctuations arc able to grow unlimitcdly and 
which will shrink (see Fig. 3). Typically, fluctuations of 
Sp H /p H ~ 0.001-0.1 above the bulk point are needed for 
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the nuclcation of drops with R ~ 10 — 1000 fm. However, 
the nucleation rates Y vary over several orders of magni- 
tude. Our results show that drops with R ~ 2 — 800 fm 
have a huge nucleation rate while those with R > 800 fm 
are strongly suppressed (see Fig. 4). 

We have also shown that fluctuations in the tempera- 
ture and in the chemical potential of trapped neutrinos 
arc not very important for dcconfincmcnt. This is in 
contrast with previous results found within the frame of 
the MIT Bag model (e.g. in [26[ it is argued that nu- 
cleation is suppressed at T < 2 MeV and in [2(| it is 
found that neutrino trapping precludes deconfinement). 
Instead, fluctuations in the energy density are the more 
efficient way to trigger the transition. 

Notice that the nuclcation rate and the typical radii 
of deconfined drops are also very different from the val- 
ues found within the MIT Bag model (see e.g. and 
references therein). The drops studied in Q have typi- 



cally radii less than 10 fm and a long nucleation time, in 
contrast ours may have much larger radii and nucleate 
almost instantaneously. This is due to the use of dif- 
ferent equations of state for the quark phase as well as 
for the different treatments of the surface and curvature 
terms. While in 0] the surface tension is assumed to be 
constant (a = 30 MeV fm" 2 ), in our work a and 7 are 
calculated self consistently within the MRE formalism 
resulting non-constant values around 140 MeV fm~ 2 and 
110 MeV fm -1 respectively. Since our surface tension 
is larger, larger critical drops are obtained. The typical 
\5Q C \ is also larger and results in huge nucleation rates. 
In the context of protoneutron stars the main conclusion 
is that if the bulk transition point is attained near the 
star centre, quark matter drops with R < 800 fm will nu- 
cleate instantaneously. Since the bulk transition density 
is ~ 5 — 6/50, this should happen for stars with masses 
larger than ~ 1.5 — 1.6M Q . 
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